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Abstract 



We make an attempt to resolve the discrepancy of the observed e + e~ mass 
spectrum in the K + — > ir + e + e~ decay with that predicted by meson domi- 
nance. To this end we investigate the properties of the p° propagator. We use 
dispersion relations to evaluate the running mass m of the p° resonance 
without adjustable parameters. To improve the convergence of the disper- 
sion integral, the momentum dependence of strong vertices is taken from the 
flux-tube-breaking model of Kokoski and Isgur. The obtained behavior of 
fripit) at small momentum squared t makes the K + — ► ir + e + e~ form factor 
rise faster with increasing t than in the original meson dominance calculation 
and more in agreement with the published data. As a consequence, the me- 
son dominance prediction of the p + p~ /e + e~ branching ratio changes slightly, 
from 0.224 to 0.236. We do not see any possibility to accommodate into the 
meson dominance approach an even steeper e + e~ spectrum, indicated by the 
preliminary data of the E865 collaboration at BNL AGS. 
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The decays K + — > n + £ + £~ (£ = e , //) have been the subject of intensive theoretical 
studies since the late fifties (see Jl]-[3| and references therein). The picture of later theoretical 
development can be gained by inspecting [|H7]] and papers cited there. The decay K + — > 
7i + £ + £~ was experimentally observed in 1975 in its e + e~ mode || and in 1997 in the fi + fi~ 
mode ||. Other experiments include a more precise measurement of the e + e~ mode by the 
BNL-E777 collaboration f!0 |, unpublished e + e~ data of the BNL-E851 collaboration |TT 



and the current BNL-E865 experiment, capable of measuring both modes (T2| with high 
precision and statistics. 

Today, it is customary to interpret the experimental results in the framework of the 
chiral perturbation theory [p|,P^,[7|. Unfortunately, this theoretical framework contains free 
parameters, one in the p A order J5], two in the p 6 order J7|. This, on one hand, diminishes the 
predictive power of the theory but, on the other hand, gives more room to experimentalists 
when trying to fit theoretical formulas. 

On the contrary, as we have shown recently ||14|| , the meson dominance offers a parameter- 



free description of the K + — > ty + £ + £~ decays. The relevant Feynman diagram is shown in 
Fig. |I|. The corresponding formula for the differential decay rate in dilepton mass M has 
the form generally expected for the one-photon approximation, namely 

dT(K+ -> 7v + £ + £~) „ „, x ,/2, 2 2 ,,2x / 4m?7 2 ™A r,,,^ 2 

-^ = C MX 3 / 2 (m 2 K+} ml +} M 2 ) \ 1 - 1 + ^ f F(M ) (1) 



dM 

with X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz and the form factor given by 

2 

777, 

F(t) = • (2) 



m p 



The normalization constant C is not given by first principles, but can be determined using 
the data other than those on the K + — > 7r + £ + £~ decays themselves. Concretely, by the 
experimental information about the r~ — > 7t~7t~tt + i> t and K + — > /i + z/ At decay. In Ref. M|, we 



used the decay rate of the r~ — >■ a^z/ T decay, the ai(1270) decay width, and the K + — > /i + z/^ 
branching fraction. In this way we obtained B(K + — > 7r + e + e _ ) w 3.1 x 10 -7 , not in 
contradiction with experiment (2.74 ±0.23) x 10~ 7 |[L5|| . The approximative character of our 
result was caused by the badly known a x decay width. 

Formula ([J) makes a definite prediction for the /i + /i _ /e + e~ branching ratio even if C is 
badly known. The number is 0.224 with an error which is negligible under the circumstances 
because the ratio is a function of the masses of participating particles only. We use it and 
the experimental e + e~ branching fraction to predict 

B(K + -> 7r+/iV~) = (6.2 ± 0.5) x 10" 8 , (3) 

in agreement with the later measurement || of (5.0 ± 1.0) x 10 -8 .[] 



1 We follow the convention of (Tf|, where the statistical, systematic and theoretical errors given in 
||] are summed quadratically. 
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The t-dependence of the form factor (0) can be, for the purpose of comparing with data, 
characterized by the slope variable 



am- »4^, (4) 

which is equal to 0.033 at t — and reaches 0.053 at the upper kinematic boundary t ~ 
0.125 GeV 2 . In Ref. |H| the data were fit by a linear approximation to the form factor 
[A = A(0)] 

F{t)=l + X-^- (5) 

with A = 0.105 ± 0.035 (stat.) ± 0.015 (syst.). 

This result became a little surprising after the experimental value of B(K + — > n + p, + pr) 
was published @. In fact, if one assumes the p/e universality and validity of fl5|), the 
values of B(K + — > 7r + e + e~), B(K + — > 7T + p + p~), and A must match together. And they do 
not match very well. Even for A = 0.055 (mean value minus both errors) the "predicted" 
value of the p + p~ /e + e~ branching ratio is equal to 0.235, which should be compared to the 
experimental 0.18 ± 0.04 (using [IS|, errors summed quadratically). The disagreement rises 



with A. The preliminary data of the E865 experiment Jl2[ indicates that the fault is not on 
the side of A. 

We therefore take for granted that the experimental value of A indicates, despite its large 
errors, that the meson dominance form factor (g) is too flat. In the following, we will try 
to find the possible origin of this discrepancy and a way to improve the situation without 
introducing unnatural assumptions and free parameters. 

When writing (fj) with form factor given by (|2|), the essential assumption was that the 
p° propagator in Fig. [I] can be written in a free- vector-particle form 

-<T + q^/ml 
t — m 2 p + xe 

where m p is the mass of the p° resonance, as it is seen in the hadronic production exper- 
iments. The general expression for the interacting-vector-resonance propagator is a little 



more complicated. It sounds (see, e.g., ||T6|) 



W t-ml(t)+im p r(t) ' 1 } 



where Y{t) is the total width of the p-resonance with off-shell mass ^/t, normalized at t = m 2 p 
to the nominal width V p . Furthermore, m 2 p (t) is the running mass squared and uj(s) is a 
complex function which reflects the properties of the one-particle-irreducible bubble. 

The propagator (||), which is usually used in meson dominance calculations, differs from 
(0) in three respects: 

1. A simplified structure of the q^q y term. This is not important, because this term does 
not contribute anyhow due to the transverse pirn vertex. 



3 



2. The absence of a finite imaginary part, which is justified, since most of our t-region 
lies below the 7T7T threshold. In a small window between the latter and the end of the 
t-interval it is negligible. Nevertheless, we will include it in what follows. 

3. The only real difference is in replacing the running mass m p (t) with the nominal mass 
m p , what is generally allowed only in a close vicinity of the resonance point. 

We will concentrate our effort on the last issue and study the consequences of replacing 
nominal mass of the p° resonance by its running mass in the denominator in Eq. (fj). To be 
more concrete, we will write the modified form factor in the form 

mj(0) 

F(t) = m 2 p {t)-t-im p T{t) ' (8) 

It follows from the causality of the propagator (|7|) that above the tttt threshold to, the 
m 2 {t) and m p Y{t) are boundary values of the real and imaginary parts, respectively, of a 
function analytic in the cut t-plane. We can therefore write a once subtracted dispersion 
relation |16| 



9/ n 9 t - m l^r°° rn o r(t') 

where symbol V denotes the principal value. To proceed further, we must find all important 
contributions to the variable width T(t). Without any doubt, we start with p° — > tt + tt~. 
Other candidates are, ordered according to rising thresholds: p° — > 7]7t + tt~ , p° — > ujtc , and 
p° — > K + K~&zK°K°. The relative importance of those channels can further be assessed 
by comparing the abundance of their isotopic companions in the r~ decays. This suggests 
that of those three, the ujtt final state will be most important, while the 7]tt + tt~ one least 
important. The results of actual calculations confirm this estimate. Furthermore, the in- 
spection of the r— decay fractions shows that there is no other important hadronic channel 
with quantum numbers of the p° meson. In addition, we assume that possible channels with 
the thresholds above the r~ mass may be neglected. The results obtained below seem to 
validate this assumption. 

Now we going to describe our calculation in more detail. Let's start with the most 
important contribution to V(t), which is the p° — > 7t + tt~ decay. We write the pirn vertex in 
the form 

v = f P uf 2 ) - p:~) , (io) 

where p* is the pion momentum in the p rest frame. Instead of the usual coupling constant 
we have introduced the strong form factor. Its momentum dependence was taken from the 



flux-tube-breaking model of Kokoski and Isgur JT7j. We thus write 



fpnAP* 2 ) = 9?™ exp | --^m f ( n ) 



with (3 = 0.4. We must confess that our original motivation for borrowing (|TlD from [17 



was technical: we just wanted to ensure a good convergence of the dispersion integral. But 
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it appeared later that the very reasonable result for m 2 p (t), which we will present below, 
could not be achieved without assuming: (|TlD or with a very different value of parameter 



p. Our opinion is now that the flux-tube-breaking model Ansatz ( ]TT| ) reflects correctly the 
real dynamics of the pinr vertex. We will use the same parametrization for all strong form 
factors. 

Using (|ToD and (JTT]) we easily arrive at the formula 



r p o_ 



■(*) 



9 piTIT V 

6/T t 



exp 



P 



*2 



6/3 2 



(12) 



where p* = — m 2 . The coupling constant was determined from the condition 



r p0 _ 



m. 



'150.7 ± 1.1) MeV 



(13) 



with the result = 41.7 ± 0.3. Formula (0) can be used, with obvious modifications 
also for p° — > KK. Here, the coupling constant can be determined from the r _ — > K~ K°v T 
branching fraction. We refer the reader to Ref . |14| . Taking into account the modifications 
connected with present usage of the momentum dependent strong form factors and assum- 
ing that p~~ and p° decay to their corresponding KK systems with the same rate, we get 
9 2 p ok+k- + 9 2 p o k oko = 28.2 ± 5.1. 

The pun vertex is taken in the formf] 



ypu 



fpUJIT iP 



e^Pp^p 



(14) 



with the same momentum dependence of the strong form factor as in ([II]). The coupling 
constants can be determined from decay rate T(uj — > 7r°7) = (7.2±0.4) x 10 -4 GeV assuming 
usual vector-meson-dominance form of coupling between p° and 7. The result is g 2 
155 ± 8. The contribution to T(t) is given by the formula 



1 pLUTY 



o(t) 



12tt t 



exp 



P 



6f3 2 



(15) 



with p* 2 = X(t,ml,mlo)/(At). 

The last contribution to T(t) we consider is the decay p° — > 7]7t + tt~ . We will considered 
as a two-step process: p° — > r]p° followed by the decay p° — > 7r + 7r~ . The mass squared of the 
parent p° is t, that of the daughter p° is s < t. Thanks to the daugher decay matrix element 
being transverse, the decay rate of the whole process factorizes into two parts |I8]|. The 



first of them is given by formula (|15|) with obvious modifications; the second one contains a 
Breit-Wigner term with decay rate of p° — > 7r + 7r~. The only new element is the prjp coupling 
constant, which is determined from T(p° — > 777) = (3.6 ± 1.3) x 10 -5 GeV as g 2 „ n = 55 ± 21. 
The contribution to T(t) is given by 



iprip 



2 The coupling constants here were made dimensionless, contrary to [14], by introducing the p 
mass in the denominator. 
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r ^-W = "36^T V + (s-m^ + m^is^X 6^-/ ^ ' (16) 



where 



A 1/2 (t, g ,m^) 
2^ 



p; = :t v . (17) 



It seems to be a sort of conundrum that the right hand side contains the same quantity, 
the contribution to which we aim to determine, namely T(s). Under different circumstances 
we would be forced to repeat whole procedure several times in search for a selfconsistent 
solution. Fortunately, here it shows that the result depends only little on the form of T(s). 
We compared the case of fixed width T p with the case of T pnw (s) and found only tiny 
differences. We picked the result of the latter choice. 

Now we have collected all pieces and can add them to form the total T(t) and evaluate 
the dispersion integral. To be sure that we have things under control, we proceeded in a less 
straightforward way. We first took the basic (p°7r + 7r~) contribution alone and determined 
m 2 (t). Then we did the same thing for the basic contribution combined with three other 
contributions taken individually and compared the changes against the basic contribution 
alone. In this way we determined the following sequence of contributions (most important 
first): cg>7t°, KK, ^7r + 7r _ . Then we started again and added the contributions cummulatively, 
in the order just shown. The resulting m^(i)'s are depicted in Fig. |[ We can see that the 
procedure of adding contributions converges to a very reasonable result: a wide plateau, 
the derivative at t — m 2 almost vanishing. Our input parameters, coupling constants, have 
relatively large errors. So it would be possible to vary them within limits in effort to find even 
better solution (characterized by the vanishing derivative at t — m 2 . Another possibility 
would be to fiddle with uj a little around the breaking-flux-tube model |T7] preferred solution 



(3 = 0.4. We made only one try in that direction. We found that the derivative vanished 
if the coupling constants squared of all three additional contributions were diminished by 
8%. But the behaviour of the running mass squared in the region which interests us most 
[0 < t < (rriK+ — rn n +) 2 } did not change by that move at all. We therefore believe that our 
determination of m 2 p (t) at low t is stable and trustable. 

Before we draw conclusions about the K + — > tt + £ + £~ form factor, we must mention 
one correction we should make in order to be consistent with the formalism we used in 
our dispersion relation evaluation of the p° running mass. We should include the same 
momentum dependence of strong form factors also into our basic diagram, Fig. [TJ. Here, it 
applies to the aipn vertex and leads to the following modification of the form factor (|8|): 



m 2 (t) — t — im p T(t) [ 48m 



K+ 



Anyhow, to see the effect of the running mass alone, in Fig. |3| we present three curves: the 
old meson dominance form factor calculated from Eq. (0) (dashed curve); the form factor 
coming from the running mass with the vertex correction ignored, Eq. (0) (dash-dotted 
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curve); and the form factor reflecting both effects, Eq. ([19]) (solid curve). The latter is what 
we consider the final product of our study. 

Fig. [| brings the same information but in a form which is better suited for comparison 
with the experimental mass spectra. It shows the dependence of the form factor squared on 
the dilepton mass. 

The form factor calculated from (|19D has a much steeper t dependence than the original 
form factor (0). It is characterized by A = 0.043 at t — and A = 0.073 at the largest t. The 
e + e~ /fi + fi~ branching ratio calculated using ([]]) with (|19|) is 0.236. Using the experimental 
branching fraction of the e + e~ mode []T5| we get a new prediction for the n + fi~ mode 

B(K+ -> 7T V>") = (6.5 ± 0.6) x 10~ 8 , (20) 

which differs only little from the original one @. The "effective A" of our form factor, 
defined as the value of A in linear parametrization (H) that leads to the same fi + fi~/e + e~ 



branching ratio, is 0.057 (in original version of the meson dominance calculation |14| it was 
0.039). From the above we conclude that the meson dominance model with our new form 
factor is consistent with the shape of the e + e~ mass distribution as measured in experiment 
0- 



A much different story is the comparison with the preliminary data [0 of the E865 ex- 
periment at the Brookhaven National Laboratory Alternating Gradient Synchrotron. Their 
10,000 K + — > 7r + e + e _ events yielded a preliminary result of the form factor parameter 
A = 0.20±0.02. If this value is confirmed, the meson dominance model of the K + — > n + £ + £~ 
decays will be ruled out, despite its success with a parameter-free calculation of the branch- 
ing fractions. 

Let us conclude with a general comment. It is a little unfortunate that the role of the 
K + — > 7t + £ + £~ decays is sometimes shrunk to a testing ground of the chiral perturbation 
theory and other clues are not followed. Here I mean mainly the importance of these decays 
in studying the behavior of the electromagnetic form factor induced by the p° resonance 
at small t. In my opinion, there are only two kinds of processes that are able to perform 
this task. Besides the K + — > n + £ + £~ decays these are the uj —>■ n°£ + £~ Dalitz decays.0 



Concerning the latter, the only e + e~ experiment performed flljj had low statistics and was 
unable to bring the mass spectrum. The fi + fi~ experiment |20[ showed that the dimuon 
mass spectrum disagreed with the vector meson dominance hypothesis. The parallel with 
K + — > n + £ + £~ is interesting. But the kaon decays we consider here are unique in populating 
mainly the region below the tctc threshold, whereas the dilepton mass spectrum of the u 
Dalitz decays spans to much higher values. The K + — > n + £ + £r decays can serve as a unique 
magnifying glass for studying behavior of the p-induced electromagnetic form factor at small 
t. 



3 The Dalitz decays of cf>, which can, in principle, serve for the same purposes, have not been 
observed yet. 
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FIGURES 




FIG. 1. Matrix element of the decay K + — > ir + e + e in the meson dominance approach. 
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FIG. 2. Running mass squared of the p meson for different inputs to the dispersion relation: 
(i) p° — > 7r + 7r _ only (dash-dotted curve); (ii) p° — ► lott added (dashed); (hi) also p° — > i^K added 
(dotted curve); (iv) the final curve (solid) after the p° — »■ i]ir + ir~ contribution has been added. 
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FIG. 3. K + — ► 7r + e + e form factor as a function of t: (i) constant p° mass, Eq. (|2|) (dashed); 
(ii) running p° mass, Eq. (Q) (dash-dotted); (in) running p° mass and the aip-K vertex correction, 
Eq. (ID (solid). 
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FIG. 4. Same as Fig. ||, but the K + — ► 7r + e + e form factor squared as a function of the 
dielectron mass. 



13 



